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Valence quark contributions for the 7iV — > Pn(1440) form factors 
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Using a covariant spectator quark model we estimate valence quark contributions to the Fj(Q 2 ) 
and f2*(Q 2 ) transition form factors for the 7./V — > Pn(1440) reaction. The Roper resonance, 
Pn(1440), is assumed to be the first radial excitation of the nucleon. The present model requires no 
extra parameters except for those already fixed by previous studies of the nucleon. Our results are 
consistent with the experimental data in the high Q 2 region, and those from lattice QCD. We also 
estimate the meson cloud contributions, focusing on the low Q 2 region, where they are expected to 
be dominant. 



I. INTRODUCTION 

The structure of the Pn(1440) resonance, which will 
be referred to as "the Roper" in this article, has been 
a long-standing mystery. Its mass and large width are 
difficult to understood in a framework of quark mod- 
els [J-Q. The decay branches imply that it is a mix- 
ture of ttN and ttttN states (dominated by the mixture 
of 7rA and aN channels). Recently, the precision data 
from CLAS at JLab [|[ (single and double pion electro- 
production), and a global analysis from MAID [(| @] in 
the region Q 2 < 6 GeV 2 have become available. The 
extracted 7iV — > Pn(1440) hclicity amplitudes can give 
hints on the internal structure of the Roper. Indeed, the 
amplitudes (and the transition form factors) extracted 
from the experiments, turned out to show a surprising 
Q 2 dependence. 

In a simple quark model picture, the Roper can be 
regarded as the first radial excitation of the nucleon. Re- 
cent experimental data at large Q 2 [f| Q support this 
picture 0, i, i, Q. In the high Q 2 region (Q 2 > 2 
GeV 2 ) the data is consistent with the predicted Q 2 de- 
pendence of relativistic constituent and light-front quark 
models 0, 0> Hi- These models, however, fail to ex- 
plain the lower Q 2 data, particularly at the photon point, 
Q = 0. This discrepancy has been interpreted as a man- 
ifestation of the missing meson cloud effects in the mod- 
els 0, [8l4Tfl| . Unfortunately, the magnitude of the me- 
son cloud contributions is less constrained, and they are 
more model dependent than the valence quark contribu- 
tions, which are dominant at larger Q 2 . (In general, the 
meson cloud contributions are expected to fall as Q 2 in- 
creases.) In this study, we estimate the meson cloud con- 
tributions using a covariant spectator constituent quark 
model which is consistent with the large Q 2 region data 
of the CLAS. The present model requires no extra pa- 
rameters except for those already fixed by the previous 
studies for the nucleon [TTj j . 

The 7 A" — > Pn(1440) transition amplitudes have been 
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studied in several ways: treating the Roper as a 
gluon system [l2j], a pure valence quark system [4. 
18 1 , and a valence quark system dressed by qq pairs 
231 ]. Relativity has proven to be very important, not 
only for a large Q 2 , but also for the region Q 2 ~ 0. 
It may even decide the sign of the amplitudes near 

Q 2 = 0, H [H 0. The nucleon t0 R °P er transi- 
tion amplitudes have been also studied using dynami- 
cal baryon-meson coupled-channel models. They usually 
treat the baryons and mesons as effective degrees of free- 
dom, and the baryons are dressed with a meson cloud 
nonperturbatively. (For details, see e.g., Refs. [13, HJ- 
[HJ and a review (H)].) The nonvalence quark degrees of 
freedom was also introduced by coupling the pion and 
other meson fields to the valence quarks |34T - [3a ]. The 
77V — > Pn(1440) transition was also studied in lattice 
QCD [13. 

In this study we compute the 7 A — > Pn(1440) tran- 
sition form factors using a valence quark model based 
on the covariant spectator formalism [38j . Relativity is 
implemented consistently. In this model a baryon is de- 
scribed as a quark-diquark system, where the diquark 
has all possible polarizations (spin state or 1) and acts 
as a spectator. The isolated quark in the baryon inter- 
acts with the photon in the impulse approximation. The 
model has been applied successfully for the studies of 
spin 1/2 and spin 3/2 low-lying baryons [lH |39M4^ ■ We 
assume that the Roper is the first radial excitation of 
the nucleon, which also ensures that it is orthogonal to 
the nucleon (valence quark) state. Using a model with 
no extra new adjustable parameters, we compute the Q 2 
dependence of the transition form factors. We find an 
excellent agreement with the data d H in the high Q 2 
region. Furthermore, we extend the model to a lattice 
QCD regime and compare with the lattice data for a 
large pion mass (m ff = 720 MeV), and find also a good 
agreement, particularly in the low Q 2 region, Q 2 < 1 
GeV 2 . Note that under the two conditions, high Q 2 and 
lattice simulations with heavy pions, meson cloud effects 
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are expected to become small or negligible [45| . The two 
agreements mentioned above, support that the present 
model can describe well the valence quark contributions. 
Encouraged by the successful features of the model, we 
estimate the meson cloud contributions focusing on the 
low Q 2 region. 

This article is organized as follows: In Sec. [TT| gen- 
eral relations among the current, helicity amplitudes, and 
transition form factors are given. Formalisms, wave func- 
tions, and explicit expressions for the transition form fac- 
tors are presented in Sec. IIIII In Sec. IIVI numerical re- 
sults are presented, and a comparison is made with the 
lattice simulation data. Furthermore, meson cloud con- 
tributions are estimated in the present approach. Finally, 
discussions and conclusions are given in Sec. [V] 



II. HELICITY AMPLITUDES AND FORM 
FACTORS 



The electromagnetic transition between a nucleon 
(mass M) and a spin 1/2 positive parity nucleon reso- 
nance N* (mass Mr) can be written using the Dirac F£ 
and Pauli F 2 type form factors, defined by the current 



J" =u R (P+) 



Mr + M' 



F* u(P. 



(1) 



where, ur (u) is the N* (nucleon) Dirac spinor, P + (P_) 
is the final (initial) momentum and q = P+ — P— . Spin 
projection indices are suppressed for simplicity. 

Usually, experimental data measured for hadron elec- 
tromagnetic structure are reported in terms of the he- 
licity amplitudes in a particular frame. The most popu- 
lar choice is the rest frame of the final state, projecting 
the current on the photon polarization states, e^, where 
A = 0, ± is the photon spin projection. In the N* rest 
frame,_the helicity amplitudes, A\i<i and S1/2, are given 



by @,iai: 



Ai, 2 (Q 2 ) = K,-{N* 1+ \\e + ■ J\N, -i), 



(2) 



Si/ 2 (Q 2 ) = IC-(N*, +\\e ■ J\N, (3) 



Here, e = \fAna is the magnitude of the electron charge 
with a ~ 1/137, and 



K 



2na 



(4) 



Wlth K = -rmr 

N* rest frame, 



|q| is the photon momentum in the 



Q 2 +Q- 
2Mr 



(5) 



where Q\ = (Mr ± M) 2 + Q 2 , with Q 2 = -q 2 . 

The helicity amplitudes, Ai/ 2 and 6*1/2, can be related 
with the form factors F? and F 2 * via Eqs. Q-© 0,0,9: 



a 1/2 (q 2 ) = tz{f*(q 2 ) + f;(q 2 )}, 

2 , K t _ t M R + M 



(6) 



s 1/2 (Q 2 ) 

where r = 



V2'"' Q 2 
-, and 



{F?(Q 2 )-tF*(Q 2 )}, (7) 



(M R +M) 2 ■ 



n = 



MbMK 



(8) 



Note that the amplitude S\/ 2 is determined by virtual 
photons and not specified at Q 2 = 0, but obtained only 
in the limit Q 2 — > 0. In this case, one has 



A 1/2 (0)=ftF 2 *(0) 
K 
75 



(9) 



S 1/2 (0) = —(M R + M)K\^(0) 



dF^ 
dQ 2 



F 2 *(0) 



(Mr + Mf 



(10) 



The inverse relations for F*(i — 1,2) in terms of the 
helicity amplitudes are: 

FKQ 2 ) = [a 1/2 (Q 2 ) + V2^±^S 1/2 (Q 2 ) 



(11) 



f;(Q 2 ) = ^YT7 { A v 2 (Q 2 ) - v^^tS 1/2 (q 2 )} 



(12) 



If the amplitudes A 1 / 2 and S%/2 are finite at Q 2 = 0, 
Eq. (fTT|) implies that 



F*(0) = 0. 



(13) 



As for F 2 *(0), there is no particular condition. 

For the Roper a fit to the data suggests a small value 
for 51/2(0) 7]. Thus, at low Q 2 , F? and F 2 * are deter- 
mined essentially by Ai/ 2 (Q 2 ). 



III. FORM FACTORS IN A VALENCE QUARK 
MODEL 

So far, a covariant spectator quark model has been de- 
veloped and applied successfully to the spin 1/2 [III |44| 
and 3/2 39-43] ground states with no radial excitations. 
In the model, a three-quark baryon is described as an iso- 
lated quark which interacts with a photon, and a specta- 
tor diquark. The wave function of a baryon is represented 
in terms of the flavor and spin states of the quark and 
diquark combined with the relative angular momentum. 
In the spirit of the covariant spectator theory the quark 
pair degrees of freedom are integrated out to form an 
on-shcll diquark with a certain mass m o . 
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A. Baryon wave functions 

To describe the momentum distribution of the quark- 
diquark system in a baryon B, we introduce a scalar 
wave function -0 b > which depends on the relative angular 
momentum and the radial excitation of the system. As 
the baryon and the diquark are on-shell the scalar wave 
function ipB can be written as a function of (P — k) 2 
where P is the baryon total momentum and k the 
diquark momentum. We represent that dependence in 
term of the dimensionless variable , 



X.B 



{M B - m D ) 2 - (P - kf 
M B m D 



(14) 



where B is the baryon index (B = N,A,N*, etc) and 
Mb the mass. 

The nucleon wave function, ^m (P, k), can be written 
in the simplest S-wave model fll| : 



^ N {P,k) 



V2 



+ (15) 



where 0°'^ represents isospin (J) or spin (S) states cor- 
responding to the total magnitude of either or 1 in the 
diquark configuration [ll]. (See Appendix lAl for the de- 
tail.) The wave function represented by Eq. ([15]). satisfies 
the Dirac equation The scalar wave function V>jv 

(S state) is given by 



m D (fa + Xn){Pi + Xn) ' 



(16) 



where xn is obtained by inserting Mb = M in Eq. (fT4j) . 
and A^o the normalization constant. In a parameteriza- 
tion where fa > fa, fa is associated with the long range 
physics, while fa the short range physics. 

In the present approach, the spin and isospin of the 
Roper are the same as those of the nucleon. Assuming 
the Roper to be the first radial excitation of the nucleon, 
we can write the Roper scalar wave function ipR as 



rl> R (P,k) = N h 



fa ~ XR 

{fa + xr) 



m D (fa +Xr)(P2 +XR.)'' 



(17) 



where fa is a new parameter which will be discussed 
later, and Nr the normalization constant. The sign of 
the normalization constant can be fixed by the exper- 
imental data Q. This particular form with the order 
one polynomial in xr ( as fa ~ Xr), is motivated by the 
harmonic-oscillator potential model for the three- quar k 
system. A similar dependence was adopted in Ref. (171 ]. 

Note that the wave function represented by Eq. (JT7J) 
preserves the short range behavior as presented in the 
nucleon scalar wave function, and that it simultaneously 
modifies the long range properties through the factor 
(ffi+xn) • ^ e R°P er wave function ^ r also satisfies the 
Dirac equation with mass Mr. 



The parameter fa in the Roper scalar wave function 
in Eq. (fl"7|). will be fixed by the orthogonality condition 
between the nucleon and the Roper states. Nr will be 
fixed by J k | , 0b| 2 = 1 at Q 2 = 0, the same as that for 
the nucleon [ll(. Thus, the description of the Roper re- 
quires no extra parameters to be fixed for the present 
purpose, since the parameters fa and fa have already 
been fixed [ill. 



B. Transition form factors 

In the covariant spectator quark model the transition 
current for the — > Pn(1440) reaction can be written 
in a relativistic impulse approximation (ill . l39l . |40| , 



\ J k 



(18) 



A 

where j k = j g^rc^p with En the diquark on-shell en- 
ergy, P+- P- = q (Q 2 
current parameterized as 



-q ), and yi is the quark 



^=ji(Q 2 b^ + j 2 (Q 2 )^. (19) 

The Dirac and Pauli (j'2) quark form factors in the 
above are also decomposed as 

h{Q 2 ) = \h+{Q 2 ) + \h-{Q 2 )r3, (< = 1,2). (20) 

The quark form factors fi± are normalized as /i±(0) = 1 
and /2±(0) = k± (isoscalar and isovector quark anoma- 
lous moments). Their explicit expressions can be found 
in Refs. [UlSliil and in Appendix E] 

Using the expressions for the nucleon and Roper wave 
functions, Eqs. (fT5 |) - ([l~7l) . and the hadronic current, 
Eq. (TT8"|) . one gets explicit expressions for F^(Q 2 ) and 
F*{Q 2 ): 

FAQ ) = 2 3lI+ 2 (M R + M) 2 + Q 2 nI 



Mr + M 



Q 2 



M (M R + M) 2 + Q 2341, 



(21) 



p;m = 7 



3 Mr + M 



. T (Mr + M) 2 . 



4 M {M R + M) 2 + 

Mr + M (Mr + Mf - 3Q 2 



+ 



2M 



(Mr + M) 2 + Q 2 



Hi, (22) 



where I(Q 2 ) is the overlap integral for the Roper and 
nucleon scalar wave functions: 



and 



J(i+2) 



1(Q 2 ) = J MP+,k)iP N (P-,k), (23) 



\u+(Q 2 )~\h-(Q 2 )Tz, (i = l,2). (24) 
6 6 
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FIG. 1: Valence quark contributions (solid lines) calculated in the present model for F± and F£. The CLAS data [j| (squares with 
error bars) and the MAID fit (dashed lines) are also shown. The result for FjT (0) is obtained from the PDG result [4(j | by converting the 
A 1/a (0). [See Eq. (5|.) 



C. Orthogonality condition between the nucleon 
and the Roper states 



CLAS data 
Spectator (valence) 
MAID 




Q 2 (GeV 2 ) 



FIG. 2: Valence quark contributions for F£, but partly magnified 
from Fig. [T] See also the caption of Fig. [T] 



Equations (f21 j) - (f22j) are the main expressions of the 
present model for the transition form factors. At Q 2 = 0, 
one has 



*?(o) = -0i+j 3 )i(o), 

= 1(1 + 7-3)1(0). 



(25) 



The orthogonality condition between the nucleon and 
the Roper states in the present approach is: 



1(0) = / i> R (P+,k)i> N (P-,k) = 0. (26) 

Jk 

This can be regarded as a generalization of the nonrel- 
ativistic orthogonality condition when Q 2 = in the fi- 
nal Roper rest frame, namely, P+ — (Mr, 0,0,0) and 



V 2M K 



,0,0,- 



2M R 



. In Appendix [B] we dis- 



cuss more on Eq. (|26l) and the nonrelativistic orthog- 
onality condition. Note that the condition Eq. (|26|) is 
not automatically verified for the nucleon and the Roper 
wave functions in Eqs. (fT6 |) -([T7 | . The orthogonality con- 
dition is satisfied only for a particular value of /3s that 
sets 1(0) = 0. 

A direct consequence of the orthogonality condition 
Eq. ([261), is that ^T(O) = 0, and also that P 2 *(0) = 0. 
F*(0) — is consistent with the general properties of 
the 7^ -t Pn(1440) reaction. F 2 *(0) = is a predic- 
tion of our model as a consequence of the orthogonality 
condition Eq. (1261) , but it also corresponds to an approx- 
imation, since the sea quark contributions are ignored in 
the present approach. A more accurate treatment would 
give a small value for (0)|. 



IV. RESULTS 



In this case, the desired result, -Fi*(0) = 0, is ensured only 
if 1(0) = 0. But this is indeed fulfilled by the orthogo- 
nality condition between the nucleon and the Roper wave 
functions. More detail will be discussed in next section. 



The nucleon and the Roper wave functions are de- 
scribed by Eq. (|T5j) with their scalar functions, Eqs. (|T6|) - 
(|17|) . The parameters for the Roper wave function are 
determined by those of the nucleon [ll|: Pi = 0.049 and 
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FIG. 3: Valence quark contributions for the A-y 12 and S\j2 hclicity 
MAID fit are also shown 0|. See also the caption of Fig.[T] 



/3 2 = 0.717, and by the orthogonality condition Eq. (j2"o| 
to give p 3 = 0.130 and N R = 3.35. 



A. Results for Fx and _F 2 * 

Before presenting the results, we emphasize again that 
there are no new free parameters to be fixed in the model. 
Thus, for a given Q 2 , the form factors can be calculated 
using Eqs. (j2"Tj) - ([2"2"j) . The valence quark contributions 
calculated in the present model for the Ff and F 2 * form 
factors are shown in Fig. [T] by the solid lines. In Fig. [5] 
also a part magnified from Fig. [T] for F 2 * is shown. The 
magnitudes of the present results are consistent with con- 
stituent and light- front quark models d, d, [0} ■ For con- 
venience, we present in Fig. [3] the helicity amplitudes 
calculated in the Roper rest frame, using the relations 
Eqs. ©-©. 

In Fig. [1] we also compare the valence quark contribu- 
tions with the CLAS data for Fj* and F 2 *. The CLAS data 
were extracted combining dispersion relations with a uni- 
tary isobar model analysis Q. As one can see our result 
is very close to the data in the region Q 2 > 1.5 GeV 2 
(see especially Fig. [5] for F 2 ), which supports the idea 
that the meson cloud contributions are suppressed in the 
high Q 2 region, and the assumption that the Roper is the 
first radial excitation of the nucleon. This achievement 
may be impressive since we have introduced no extra new 
parameters as already mentioned. One only may identify 
the long range behavior of the Roper wave function with 
that of the nucleon, and may ensure the orthogonality 
between the Roper and the nucleon wave functions. In 
Fig. [T]we also show a fit to the data from Ref. [6j, us- 
ing the code MAID2007 (abbreviated MAID) Q. Both 
the CLAS data and the MAID fit are similar for Q 2 > 2 
GeV 2 region, although some differences may be notice- 
able, in particular for F* . This discrepancy in F x * can be 




amplitudes in units of 10~ 3 GeV -1 / 2 . The CLAS data [| and the 



also a consequence of the different data sets used in the 
analysis between the CLAS data and the MAID fit. Our 
predictions for F* are closer to the CLAS analysis than 
to the MAID fit, although the differences are comparable 
with the error bars of the CLAS data. 



Next, we discuss the asymptotic behavior forK* and 
F 2 * as Q 2 -> oo. Perturbative QCD (pQCD) g3,li| pre- 
dicts F^(Q 2 ) ~ 1/Q 4 and F 2 *(Q 2 ) - 1/Q 6 as Q 2 -> oo, 
apart from the logQ 2 corrections. The predictions of the 
present model from Eqs. (f2~Tj) - (l2"2"l) . are consistent with 
these results: F^(Q 2 ) ~ 1.211 and F 2 *(Q 2 ) ~ 13.1^-, 

with 1 = ©(log Q 2 )/Q 4 . (See Appendix G of Ref. g| 
for details.) For the helicity amplitudes, our results are 
also consistent with the pQCD predictions 1 : A 1 / 2 ,<S , 1 / 2 ~ 
1/Q 3 [13, Hil, again apart from the logarithmic correc- 
tions. Note, however, that based on what observed for 
the nucleon [ll|, and especially for the 7 A — > A tran- 
sition form factors 40], the predicted scaling behaviors 
may not be observed in such a small Q 2 region. Also it is 
not realistic to expect that our calibration of the quark 
current based on the reactions in the regime Q 2 < 6 
GeV 2 , can be naively extended to the high Q 2 region 
where pQCD is valid (such as Q 2 - 100 - 1000 GeV 2 ). 



From Refs. [47],[4|] one has 

G+ ~ 1/Q 3 , Go ~ 1/Q 4 , (27) 

where G+ = A1/2 an d Go = t^jt<Si/2' Note that the extra power 
in Go in contrast to G_|_. This takes account of changing the 
helicity one unit between the initial and final states |47| . l48ll . 
Thus, one gets, A x / 2 ~ 1/Q 3 and 5x/ 2 ~ 1/Q 3 - 
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FIG. 4: Valence quark contributions in the lattice regime (dashed lines) corresponding to the pion mass m n = 720 MeV. The lattice data 
with = 720 MeV (diamonds) are from Ref. [37j . The solid lines are the valence quark contributions calculated with the physical pion 
mass, m w = 138 MeV. The mass values used corresponding to the lattice regime (m w = 720 MeV) for the nucleon, Roper and p meson 
are, M = 1.48 GeV, M R = 2.53 GeV, and m p = 1.083 GeV, respectively. 



B. Comparing with the lattice results 

Extending the model to the lattice simulation regime, 
we compare our results with the lattice QCD results. 
Hereafter, the results obtained in the lattice regime will 
be referred to as "lattice regime" , which will be explained 
below. One can expect that, the heavier the pion mass 
becomes in the lattice simulations, the closer our results 
become to the lattice one, since the meson cloud effects 
become smaller. A comparison is made in Fig. H]for the 
lattice data corresponding to m„ = 720 MeV [37[. To 
extend the model to the lattice regime, we include an im- 
plicit dependence on the pion mass for the hadron masses, 
following the procedure proposed in Rcfs. [4^, [H3] - This 
is done by replacing the hadron masses in the model by 
the respective lattice masses in the baryon wave func- 
tion and the quark electromagnetic current [Eq. (fT9"]) ] . In 
the lattice regime (the dashed lines in Fig. 2]) we have 
used the nucleon and the Roper masses, M = 1.48 GeV 
and Mr = 2.53 GeV, respectively, corresponding to the 
pion mass = 720 MeV |37|. We have also used the 
value, 1.083 GeV for the p meson mass, according to the 
parameterization in Ref. [5lj ]. 

One can see in Fig.@]that the lattice regime results are 
fairly consistent with the lattice data for F 2 *. As for F*, 
the lattice results overestimate our predictions for Q 2 > 
1.5 GeV 2 , as well as the covariant spectator model results 
"Spectator (valence)" [the same as those shown in Figs.Q] 
and |2]. However, for Q 2 < 0.6 GeV 2 , the agreement is 
excellent. More lattice QCD simulation data with smaller 
pion masses are desired to constrain better the valence 
quark contributions in the low Q 2 region. The agreement 
between the lattice regime results and the lattice QCD 
data, supports that the present estimate of the valence 
quark contributions is reasonable. 



C. Estimating the meson cloud effects 

To estimate the meson cloud effects, we take the va- 
lence quark contributions estimated in Sec. II V Al as a ref- 
erence. The valence quark contributions will be denoted 
by F^ are (i — 1,2). We need to know a full contribution 
to estimate the meson cloud contributions. As a first 
step approach, we regard the MAID fit results [6[ , which 
will be denoted by F t MAID {Q 2 ){i = 1,2), as the total 
contribution for the form factors. Then, the meson cloud 
contributions i 7 ™ , may be estimated by 

F™ c {Q 2 ) = F l MAID {Q 2 )-F* are {Q 2 ) 1 (t = l,2). (28) 

The results for F™ c are presented in Fig. [5] To have 
an idea for the uncertainties in this estimate we use an 
analytical form for the each form factor upper limit of 
the CLAS data error bars [denoted by a(F*)(i = 1,2) 
below], and calculate the bands with one-standard devi- 
ation. Specifically, we use <r(F*) ~ 0.015(1 - e~ 3Q ) and 
a(F 2 *) ~ 0.05- 0.009Q 2 . 

First, we analyze F*. A situation for Fj* is a bit sub- 
tle. The magnitude of the data is smaller than that of 
the F£ . The upper limit for \Fi\ is around 0.1, which 
is about half of that for [F^l- The meson cloud con- 
tributions are negative, and they can be between 0.01 to 
0.03 in the high Q 2 region. Probably, the contributions of 
0.03 are overestimated, since they are expected to become 
small in the high Q 2 region. We note however, the meson 
cloud contributions here arc estimated by subtracting the 
model results from the MAID fit. There are differences 
in the CLAS data and the MAID fit in the high Q 2 re- 
gion, since the MAID analysis [5j uses a different data set 
from that of the CLAS @. The clarification between the 
differences in the two data sets, and more precision data 
in this high Q 2 region (less contaminations from meson 
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cloud), are desired to constrain better the meson cloud 
contributions for F*. 

Next we discuss F| . Because of the larger magnitude 
of F£ , the contributions from the meson cloud may be 
easier to estimate. The meson cloud contributions for 
F 2 * are negative as well as those for F*. As one can 
see in Fig. [5] the meson cloud gives main contributions 
in the low Q 2 region. (At Q 2 = 0, they are the only 
contributions.) However, as expected, the meson cloud 
contributions become small for Q 2 > 1.5 GeV 2 , and the 
valence quark contributions dominate. In higher Q 2 re- 
gion, the magnitude of the meson cloud contributions are 
of order of the error bars (~ 0.01) of the data. 

Overall, our estimate of the meson cloud contributions 
qualitatively agree with those of the recent works [2(J [H, 
l29j , particularly for Fj* . 

In terms of the helicity amplitudes the meson cloud 
contributions at Q 2 = are estimated by A™ c 2 (0) — 

(-61.4 ±4.9) x 10~ 3 GeV~ 1/2 and S^(0) = (-40.2 ± 

2.6) x 10~ 3 GcV -1 / 2 . Here the meson cloud is the only 
contribution for A 1 / 2 (0). As for £1/2(0), the meson cloud 
contributions cancel significantly with those of the va- 
lence quark. This feature can be understood better by 
observing Fig. [3] 

Thus, we predict negative contributions from the me- 
son cloud for A 1 / 2 and S\/2 in the low Q 2 region. This 
is consistent with the estimate made in Ref. [36| . where 
meson cloud contributions are at least 70% in A 1 / 2 (0) 
(100% in our case), although in their case the bare con- 
tributions differ in sign from our results. A similar dif- 
ference exists also for Sin, in contradiction to the CLAS 
data in the high Q 2 region. Also Ref. [30] predicts signif- 
icant negative contributions from meson cloud for S1/2 
in the low Q 2 region, implying large negative values for 
S1/2 in that region, contrary to what is suggested by the 
CLAS data 0]. 



We also attempt to compare our results with those of 
the recent calculations based on the dynamical meson- 
baryon coupled-channel model of Sato and Lee [To|, 
[27l ] - There, the bare contributions are very small for 
^1/2(0) E3- In the region, Q 2 < 1 GeV 2 , the abso- 
lute magnitudes of the meson cloud contributions (their 
amplitudes are complex numbers) are dominant in both 
amplitudes, although the meson cloud contributions de- 
crease in the high Q 2 region for Ai/ 2 27]. Unfortunately, 
their analysis [28[ leads to the amplitudes different from 
those of the CLAS. Thus, it is expected that their meson 
cloud contributions and ours differ quantitatively. 



V. CONCLUSIONS 

In this work we have studied the transition form factors 
for jN — > Pn(1440) reaction using a covariant specta- 
tor quark model. The model is fully relativistic, and has 
no new adjustable parameters. This is fulfilled by the or- 
thogonality condition between the nucleon and the Roper 
states. The model can describe the reaction well in the 
high Q 2 region (1.5 GeV 2 < Q 2 < 4.0 GeV 2 ), supporting 
the idea that the Roper is the first radial excitation of the 
nucleon. In the low Q 2 region (Q 2 < 1.5 GeV 2 ) the model 
predictions deviate from the experimental data, and this 
fact suggests that the meson cloud contributions (missing 
in the model) are significant in this Q 2 region. 

To check whether or not the description of our esti- 
mate for the valence quark contributions is realistic, we 
have extended the model to the lattice regime to com- 
pare with the heavy-pion lattice QCD data. Our results 
are in an excellent agreement with them in the low Q 2 
region. This fact supports that the valence quark de- 
grees of freedom is well described in our model. Having 
a good confidence in the description of the valence quark 
sector, we have estimated the meson cloud contributions 
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for the Dirac (Pf*) and Pauli (P 2 *) type transition form 
factors. A characteristic feature of the model, F£ (0) = 0, 
is a direct consequence of the orthogonality condition be- 
tween the nucleon and the Roper states. However, this 
may not have significant consequences in the estimate 
of the meson cloud contributions, since they appear to 
be dominant at low Q 2 , or conversely, the valence quark 
contributions are expected to be small at Q 2 = [J]. 

In this work we have shown that the transition form 
factors Fi and P 2 * for the jN -» Pn(1440) reaction may 
be more suitable quantities to study than the helicity 
amplitudes in order to disentangle the meson cloud con- 
tributions. For example, \Fi\ is small for both in experi- 
mental data and the valence quark model because of the 
cancellation between the A 1 / 2 and 5i/ 2 amplitudes. This 
implies also the meson cloud contributions are small. As 
for P 2 , it is mainly determined by Ai/ 2 (Q 2 ) in the very 
low Q 2 region, which measures the meson cloud contri- 
butions in our approach (for Q 2 — » 0). 

We admit that our estimate of the meson cloud con- 
tributions has some limitations due to the parameteriza- 
tion of the amplitudes A1/2 and S1/2 based on the MAID 
analysis. Furthermore, we have faced that there are some 
discrepancies between the two data sets of the CLAS and 
MAID in the high Q 2 region. More precision experimen- 
tal data in this high Q 2 region are desired to constrain 
better the meson cloud contributions. 

For a possible improvement, there is a room to add one 
extra free parameter to the model. The parameter j3\ 
(fixed by the nucleon wave function) in the factor ^~x« 
in the Roper wave function, can be replaced by a new free 
parameter P\ to introduce a new long range scale, and 
may be adjusted by more precise, high Q 2 data. Needless 
to say, more lattice QCD data are also useful to constrain 
the valence quark contributions better. 

We plan to extend the present model for higher mass 
resonance region, like for instance, 5*11(1535). In the past 
dynamical coupled-channel models are very successful in 
the description of the meson-baryon electro- and photore- 
action, involving resonances such as A(1232), Pn(1440) 
and 5*11(1535). However, such models require a parame- 
terization in the interaction with the quark core. In gen- 
eral, it is based on the baryon-meson phenomenology and 
not based on the quark (and gluon) degrees of freedom. 
The approach developed in this work is very promising 
to study the valence quark contributions in the meson- 
baryon systems. It can be tested, or compared with the 
bare parameters used in the dynamical coupled-channel 
models. An independent test may be to compare with the 
'bare' contributions (no meson cloud) determined by dy- 
namical meson-baryon coupled-channel models by fitting 
e.g., to the cross section data at each Q 2 |26l428| . This 
kind of an extraction of the 'bare' contributions was done 
for the 7 AT — > A reaction I26j , as well as other resonances 
including the Roper [27l l28| . However, often the meson 
cloud contributions become larger than the quark core 
contributions, and an estimate of the bare contributions 
are sometimes difficult to get as a smooth function of 



Q 2 . Another possible method is to use the lattice QCD 
data in the heavy-pion regime to fix the valence quark 
contributions. This method was successfully applied in 
Ref. [10] for the nucleon to A electromagnetic transition. 
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Appendix A: Covariant spectator quark model 

Below, we present some details of the covariant spec- 
tator quark model. 



1. Wave functions 

In the covariant spectator quark model the S-state 
wave function for the nucleon with the (initial) momen- 
tum P_ is given by [111 ]. 



1 
71 



6?m(P_ 



}e P _U a (P^) i> N {P-,k), 
(Al) 



where k is the diquark momentum and (j) / 1 is the com- 
bination of the quark flavors associated with the isospin 
or 1 diquark. The spin-1 polarization vector is repre- 
sented by £p_ in a fixed-axis base [52[ , and u and U a are 
the spinors related by [39j . 

U a (P) = ± l5 (j a -^ju(P). (A2) 

The Roper wave function with the (final) momentum 
P + is defined by 



1 

V2 



$u R {P + ) - <f>Wp + U' a (P+)} 4>R(P+,k), 

(A3) 

where ur and U' a are spin states associated with a spin 
1/2 particle and mass Mr, and e'jZ is the polarization 
vector. 

In both cases the scalar functions, tpisi{P-,k) and 
■0i?(P+, k), are the functions of (P_ — k) 2 and (P+ — k) 2 , 
respectively. 
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re+ k_ c+ c_ d+ ci_ A 
1.639 1.823 4.16 1.16 -0.686 -0.686 1.21 



TABLE I: Quark current parameters. 



2. Quark form factors 

The quark current associated with Eqs. (TT9)) and (|20[) is 
expressed in terms of the quark form factors fi±{i = 1, 2), 
inspired by a vector meson dominance form: 



/i±«t) = a + (i-a) ., 

/ 2± (g 2 ) = K± (d ± 



MlQ 2 



: (M 2 + Q 2 ) 2 
+ (i-d±) 



M 2 + Q 2 



(A4) 
,(A5) 



In the above, A defines the quark charge in deep inelas- 
tic scattering, k± are the isoscalar and isovector quark 
anomalous magnetic moments. The mass m v (M/j) cor- 
responds to the light (heavy) vector meson, and c±, d± 
are the mixture coefficients. In the present model we 
set m v = m p [k, mj) for the light vectorial meson and 
Mh = 2M (twice the nucleon mass) to represent the short 
range physics. The values of the parameters we re p revi- 
ously fixed by the nucleon elastic form factors [ll[, and 
the values are presented in Table HI Note that the present 
model uses cL_ = d_ . 



3. Quark form factors and asymptotic expressions 

In the present model, the asymptotic expressions for 
Q 2 — > oo, associated with the quark current of Eqs. (|20|) 
and (|24|) are given by 



h - 






(A6) 


33 - 






(A7) 


h - 




-3*-)^, 


(A8) 


34 ^ 




^ T 

K ->Q 2 ' 


(A9) 



with 



(1 



d+)Ml 



(A10) 



where, we have used the relation, d+ = d-, correspond- 
ing to the parameterization of the model II in Ref. [j"lj |. 
In particular, T — 5.54 GeV 2 should be noted. 

Using Eqs. (|A6[) - (|A9| we can derive the expressions for 



F* and F* for Q 2 -> oo. 



Appendix B: The orthogonality of the nucleon and 
the Roper states 

In a static, nonrelativistic formalism the wave function 
is a function of the particle's three-momentum, k. In the 
limit where there is no momentum transfer, the argu- 
ments in the initial (ipi) and final (ipf) wave functions 
are the same. This corresponds to q = k — k = 0, or 
Q 2 = — q 2 = 0. Then, the overlap integral between the 
orthogonal states, X, is given by 



rf 3 k 

— W ^.(k)=o. 



(Bl) 



The equivalent relation for the covariant spectator the- 
ory for Eq. (|Blj) is 



= o, 



(B2) 



3 2 =0 



where the subindex Q 2 = indicates that Q 2 = 
— (P + — P_) 2 = 0. The simplest case is in the initial (or 
final) baryon's rest frame. For the equal mass case, 
P+ = P- = (M, 0,0,0). 

Next, consider the inelastic case with the masses of the 
initial and final states, M and Mr, respectively. The 
generalization of the condition Eq. (|B2[) would corre- 
spond to P_ = (M, 0,0,0) and P+ = (Mr, 0,0,0), with 

p+-p- = ( M *X~° 2 >°' °> iii) when i = (iii 



q 



is given by Eq. JSJ.) This gives also 



Q 2 = Q* 



-(Mr — M) 2 



(B3) 



which will be denoted by the pseudo-threshold point, the 
point where both initial and final state are at rest. This 
situation is, however, unphysical for Mr ^ M because 
Q 2 < 0. 

The generalization of Eq. (|B2|) for unequal mass case 
would be, 



r f (P+,k)MP-,k) 



= 0. 



(B4) 



Q 2 



Since the physical reaction are restricted to Q 2 > 0, we 
need to redefine the orthogonality condition in the un- 
equal mass case. To do so we impose the condition, 



r f (P+,k)MP-,k) = 0, 



(B5) 



where the four-momenta P+ and P_ are respectively de- 
fined by, 



{Ml + M 2 Ml-M 2 

P- = — ,0,0, ^ 

V 2M R 2M r 

P+ = (Mr, 0,0,0), 
which correspond to Q 2 = 0, but q 2 ^ 0. 



(B6) 
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The use of Eq. (|B5I) may be regarded as an approxima- 
tion, or the simplest relativistic extension for the orthog- 
onality condition. The exact treatment needs to impose 
the overlap integral to vanish at the pseudo-threshold 
as in Eq. (|B4I) . However, it would require an exten- 



sion of the wave functions to the unphysical region, and 
beyond a scope of the present study. Instead, we use 
Eq. (|B5|) . Both prescriptions should give similar results 
when (Mr — M) 2 is small enough. 
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